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ABSTRACT 

We investigate ^^-interactions of free nucleons and their impact on the electron frac¬ 
tion (Yg) and r-process nucleosynthesis in ejecta characteristic of binary neutron star 
mergers (BNSMs). For that we employ trajectories from a relativistic BNSM model 
to represent the density-temperature evolutions in our parametric study. In the high- 
density environment, positron captures decrease the neutron richness at the high tem¬ 
peratures predicted by the hydrodynamic simulation. Circumventing the complexities 
of modelling three-dimensional neutrino transport, (anti)neutrino captures are pa¬ 
rameterized in terms of prescribed neutrino luminosities and mean energies, guided 
by published results and assumed as constant in time. Depending sensitively on the 
adopted Ue-^e luminosity ratio, neutrino processes increase Yg to values between 0.25 
and 0.40, still allowing for a successful r-process compatible with the observed solar 
abundance distribution and a signihcant fraction of the ejecta consisting of r-process 
nuclei. If the luminosities and mean energies are relatively large compared to the 
z/g properties, the mean Yg might reach values >0.40 so that neutrino captures seri¬ 
ously compromise the success of the r-process. In this case, the r-abundances remain 
compatible with the solar distribution, but the total amount of ejected r-material is 
reduced to a few percent, because the production of iron-peak elements is favored. 
Proper neutrino physics, in particular also neutrino absorption, have to be included 
in BNSM simulations before hnal conclusions can be drawn concerning r-processing 
in this environment and concerning observational consequences like kilonovae, whose 
peak brightness and color temperature are sensitive to the composition-dependent 
opacity of the ejecta. 

Key words: nuclear reactions, nucleosynthesis, abundances - neutrinos - stars: neu¬ 
tron - dense matter - hydrodynamics 


1 INTRODUCTION 

The r-process, or rapid neutron-capture process, of stellar 
nucleosynthesis is invoked to explain the production of the 
stable (and some long-lived radioactive) neutron-rich nu¬ 
clides heavier than iron that are observed in stars of var¬ 
ious metallicities, as well as in the solar system (for a re¬ 
view, see 


Arnould et al. 2007). Despite important effort 


to model potential r-process sites, all the proposed scenar¬ 
ios face serious problems and the site(s) of the r-process is 
(are) not identified yet. Until now, type-II supernovae or 
7 -ray bursts models have failed to provide convincing ev¬ 
idence for a successful r-processing that could signihcantly 
contribute to the Galactic enrichment in r-material (Wanajo 


et al. 2011 Janka 2012 Burrows 2013). Only magneto- 


rotational supernova explosions with extremely strong pre¬ 
collapse magnetic fields and fast rotation seem to provide 
favourable conditions for r-processing, but such rare events 
are not expected to be at the origin of the global galactic en¬ 
richment in r-process nuclei ([Winteler et al.|2012 Nishimura 


et al.||2015 Wehmeyer et al.||2015). 


For this reason, special attention is now being paid to 
neutron star (NS) mergers following the conhrmation by 
hydrodynamic simulations that a signihcant amount of r- 
process enriched material, typically about 10“^ to a few 


10 ‘^Mq, can be ejected ( 

Rosswog et al.|l999 

Freiburghaus 

et al.| 199911 Arnould et al. 

2007||Metzger et al. 

2010||Roberts 


et al.| 2011| Goriely et al.||2011 

[Korobkin et al.| 

2012 

Bauswein et al.|2013| [Goriely et al. 

20 13[ [Wanajo et al. 

2014 
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2 S. Goriely et al. 


Perego et al.|[20T4 Just et ^[2015 Sekiguchi et al.||2015 |. 


Recent nucleosynthesis calculations by Just et al. (2015) 


show that the combined contribution of both the dynam¬ 
ical (prompt) ejecta expelled during binary NS or NS-black 
hole (BH) mergers and the neutrino and viscously driven 
outflows generated during the post-merger remnant evolu¬ 
tion of relic BH-torus systems can lead to the production 
of r-process elements from mass number A > 90 up to tho¬ 
rium and uranium. The corresponding abundance distribu¬ 
tion reproduces the solar distribution extremely well and 
can also account for the elemental distributions observed in 
low-metallicity stars (Roederer |2011||Roederer et al.|2012|). 
Furthermo re, recent studies ([Matteucci et al.|2014| [Komiya 


to large radii before the neutrino emission of the compact 
merger remnant becomes sizable and neutrino absorption 
can affect the electron fraction signihcantly. Absorption of 
neutrinos radiated by the massive merger remnant is also 
diminished because the tidal ejection in Newtonian models 
happens preferentially in the orbital plane (see e.g. Rosswog 
et al. 1999 Korobkin et al.|2012 ) while the neutrinos are pre¬ 
dominantly emitted perpendicular to this plane (Rosswog & 


et al.|20i4l|Mennekens Vanbeveren |2014||Stien et al.|201^ 

de Voort et al.|2014| [Vangioni et al.|2014| [Wehmeyer et] 


al. 2015) have reconsidered the galactic or cosmic chemi¬ 


cal evolution of r-process elements in different evolutionary 
contexts. Although they do not converge towards one unique 
quantitative picture, most of them arrived at the conclusion 
that double compact star mergers may be the major pro¬ 
duction sites of r-process elements. 

Despite the recent success of nucleosynthesis studies for 
NS mergers, the possibility of r-processing in these events 
is still affected by a variety of uncertainties. In particular, 
the impact of neutrino interactions is not yet studied and 
understood in detail, the main reason of which is the not 
yet manageable computational complexity associated with 
neutrino transport in a generically three-dimensional (3D), 


Liebendorfer||200^ Dessart et al.|[200^ Perego et al.||20l'4 )^ 
However, the described situation applies well only for New¬ 
tonian mergers. The situation is different in the relativistic 


Relativistic simulations of symmetric mergers do not 
show the development of prominent tidal arms, and the ejec¬ 
tion of unshocked matter is therefore not important. Instead, 
the collision shock that builds up at the interface of the 
two NSs is typically stronger than in Newtonian conditions 
(see e.g. Bauswein et al.|2013 ), leading to potentially higher 


temperatures and higher, faster rising neutrino luminosities 
( Wanajo et al.|[2QT4 Sekiguchi et al.|[2015 ). The ejecta are 


expelled fairly spherically instead of equatorially (Bauswein 
et al.|20l3 Hotokezaka et al.|2013 ) and consist of two main 


components, namely a hrst one in which very hot material is 
squeezed out from the collision interface of the two merging 
bodies, and a second, slightly delayed one that is expelled 
in waves from a torus-like belt of matter around the high- 
density core of the merger remnant. This torus is heated 
by spiral shocks that are sent out from the aspherical, wob¬ 
bling, and rotating high-density core and that also lead to 


velocities and rapid changes in time. A computationally sim- 

outward acceleration of parts of the torus matter (Bauswein 

pie and much more efficient alternative to solving the time- 

et al. 

2013 

Hotokezaka et al.|20I3||Just et al.|20I5||Wanajo 

dependent transport equation for neutrino distributions in 

et al. 

2014 

Sekiguchi et al.||20I5). 

the six-dimensional phase space is the use of a neutrino leak- 

Indeed 

, recent relativistic NS-NS merger simulations 


age scheme, in which the local neutrino net-emission (i.e., 
emission minus absorption) rate is estimated by a weighted 
interpolation between the pure emission rate and an optical- 
depth dependent diffusive loss term (e.g., Ruffert et al.|19% 
Rosswog k, Liebendorfer|[2003 ). However, neutrino absorp- 


that took into account neutrino emission by means of a leak¬ 
age scheme and absorption by an additional approximate 


transport treatment based on a moment formalism (Wanajo 
eTaTIpoTI [Sekiguchi et al.||2015| ), found that neutrino 


tion cannot be straightforwardly included in a self-consistent 
manner in a leakage scheme because information about the 
local neutrino densities is missing in such a treatment. 

It was long believed that neutrino interactions could 
not, at least not drastically, affect the initial neutron rich¬ 
ness of the ejecta. Such expectations were based on numer¬ 
ical merger models. The temperatures of the merging NSs 
in these simulations, and therefore the neutrino production 
rates, remain low until the NSs collide with each other, and 
they rise in an increasing volume only gradually on a time 
scale of several milliseconds after the hrst contact of the 


two NSs (Ruffert & Janka 1999 Rosswog & Liebendorfer 
2003 Wanajo et al.|20I4 ). At this time a signihcant fraction 


of the ejecta material is already being expelled with nearly 
relativistic velocities. Newtonian merger studies suggested 
that a large part of the ejecta material even from symmetric 
mergers (i.e., for two equal-mass or close to equal-mass NSs) 
is thrown out by tidal forces in extended spiral arms forming 
from matter of the outer faces of the merging objects during 
their hnal approach and collision ( Rosswog et al.||I999 Ko- 


teractions with free nucleons can signihcantly increase the 
electron fraction in the dynamical ejecta for cases in which 
the collapse of the merger remnant to a black hole is delayed 
or does not happen. Under such conditions nuclei with mass 
numbers A < 140 can also be created in the dynamical ejecta 
in addition to the heavy r-process elements {A > 140). Weak 
interaction processes of free nucleons can consequently affect 
the strength of the r-process and the emerging abundance 
distribution. 

The accurate inclusion of neutrino interactions in hy- 
drodynamical simulations remains a highly complex task. 
This motivated us to conduct a simple, parametric study 
in order to quantify the potential impact of weak inter¬ 
actions on the electron-tract ion evolution in merger ejecta 
and thus to explore the consequences of charged-current 
neutrino-nucleon reactions for the nucleosynthesis and possi¬ 
ble r-processing in these ejecta. More specihcally, we inves¬ 
tigate the influence of /^-interactions of electron neutrinos 
(ue) and electron antineutrinos (ue) with free n and p and 
of their inverse reactions, 


robkin et al.||20I2 ). By the tidal stretching these arms nat¬ 
urally remain unshocked and thus stay cool. Therefore they 
reach low densities so quickly that electron and positron cap¬ 
tures cannot become efficient. Moreover, the ejecta escape 


+ n ^ p + e 
+ p ^ n + 


( 1 ) 

( 2 ) 


on the Ye distribution and r-process nucleosynthesis at con¬ 
ditions representative of the dynamical ejecta expelled by 
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Weak interactions in neutron-star mergers 3 


hydro dynamical forces during NS-NS mergers. These reac¬ 
tions have been neglected in all previous studies of r-process 
nucleosynthesis for such ejecta except those of |Wanajo et al.| 
(2014), where their quantitative effects may depend on the 
adopted equation of state ( Sekiguchi et al.|[2015 ). The role 
of weak interactions for the electron fraction and the cor¬ 
responding implications for r-process nucleosynthesis, how¬ 
ever, demand further exploration in more detail, in particu¬ 
lar also by basic, parametric modeling, because a multitude 
of uncertainties will prevent rigorous, self-consistent solu¬ 
tions of the full problem in the near future. Such uncertain¬ 
ties are associated with, for example, the extreme complexi¬ 
ties of 3D energy-dependent neutrino transport in relativis¬ 
tic environments, with the neutrino opacities of dense, po¬ 
tentially highly magnetized matter, and with neutrino-flavor 
oscillations at rapidly time-variable, largely aspherical con¬ 
ditions of neutrino emission. 

The main objective of our present work is a sensitivity 
study by the use of a parametric approach. It is intended to 
motivate further explorations of neutrino effects in relativis¬ 
tic NS-NS mergers in more detail and breadth. To this end, 
we set up a simplified and idealized theoretical framework to 
test the individual roles of the different weak interaction pro¬ 
cesses, considering the density and temperature evolution of 
fluid elements ejected from a prototype hydrodynamical, rel¬ 
ativistic NS-NS merger model. For our parameter study we 
make assumptions about the neutrino emission properties 
that are guided by data taken from the literature. The elec¬ 
tron fractions resulting from the neutrino-processing of the 
ejecta elements are then used as input for nuclear network 
calculations, allowing us to immediately link the effects of 
neutrino interactions to the final heavy-element production. 

In Sect, [^the employed merger model and our treat¬ 
ment of the weak neutrino reactions with free nucleons are 
described. The effects of /3-processes on the electron fraction 
are reported in Sect.[^ and the subsequent r-process nucle¬ 
osynthesis is analysed in Sect. ^ Conclusions are drawn in 
Sect, m 


2 WEAK INTERACTIONS OF FREE 
NUCLEONS IN MERGER EJECTA 

We adopt the density evolution of ejecta fluid elements from 
a representative NS-NS merger model, namely the sym¬ 
metric 1.35M0-1.35M0 binary model obtained with the 


temperature-dependent DD2 equation of state (Hempel & 
Schaffner-Bielich |2010[|Typel et al.|2010 ) (essentially iden- 


tical to the one in |Bauswein et al.||2013| but with a higher 
resolution of particles). This relativistic hydrodynamic 
simulation also provides the temperature evolution. We do 
not apply any temperature post-processing as in |Goriely et| 
al. (2011) to disentangle temperature jumps in shocks from 
artificial heating associated with the use of a numerical vis¬ 
cosity in the smoothed-particle hydrodynamics scheme. 

For each trajectory, we follow the expansion starting at 
a fiducial density of peq = lO^^gcm”^, where we assume 
equilibrium to hold between electrons, positrons and neu¬ 
trinos for the given total lepton number provided by our 
NS-NS merger model. Below peq = lO^^gcm”^, electron, 
positron and electron neutrino and antineutrino captures 
are systematically included. Reactions of neutrinos on nu¬ 


clei are, however, neglected. As long as the temperature re¬ 
mains in excess of typically T > 10^° K, the abundance 
of heavy nuclei is determined by nuclear statistical equi¬ 
librium (NSE) at the given electron fraction, density and 
temperature. From the density peq down to the density pnet, 
at which the temperature reaches 10 GK and the full reac¬ 
tion network is initiated, the /3-interactions of free nucle¬ 
ons may affect the electron fraction Ye. If a trajectory stays 
cooler than 10 GK below the density peq, the network cal¬ 
culation is started at the neutron-drip density pdrip, i-e., 
Pnet = Pdrip — 4.2 X 10^^ gcm“^. The considered /3-reactions 
involve free nucleons, whose abundance variations are given 
by 


dt 

dY^ 


-x+Y^ + A_y;, 


= A. 


— \ 

- -I ri p ? 


( 3 ) 


dt 

where A+ = + Ag+ and A_ = Ap^ + Ag-. The A^c de¬ 

note capture rates of species x G {e“, e^, z^e, onto free 
nucleons according to the ^d-reactions, Eqs. ( |1|2[ ). The free 
neutron and proton numbers are related to Ye by 


Yl = 1-Y,-Y,NY{Z,N)., 

Z^2 

= Ye-Y,ZY{Z,N), 


( 4 ) 


where Y = A/A is the molar fraction (and X the mass 
fraction) of the nucleus {Z,N) of atomic mass A = Z N. 
Assuming that the NSE molar fractions of nuclei remain 
constant over the time step At, the time evolution of Ye can 
be related to the change of the number of free protons, and 
written as 


dYe 

dt 


= -AtotW + A* 


where Atot = A+ + A- and 


A* = A+ 


1-Y.ny 


+ X-Y,zy. 


( 5 ) 


( 6 ) 


If only free neutrons and protons are present. A* = A+. 
The impact of heavy nuclei is to increase A* towards A_ 
(for matter made of a-particles only. A* = Atot/2). The a- 


effect ( McLaughlin et ^[1996 Meyer et al. 1998| Pllumbi 


|et al. 2014), or more generally, the effect of heavy nuclei 
in binding neutrons and protons inside nuclei, is known to 
be responsible for driving Ye towards 0.5 and is included in 
this term A* . If we assume that Atot and A* remain constant 
over the time step At (or, specifically, that At is chosen such 
that Atot and A* as well as the abundance of nuclei remain 
essentially constant during the time step), Eq. © can be 
integrated analytically leading to 

Ye{t + At) ~ + V- [l - j , (7) 

Atot L -I 

where Atot and A* are estimated at time t. This equation 
is used to follow Ye from the initial density peq down to 
the density pnet, at which the temperature has dropped to 
10 GK (or, alternatively, to the drip density if temperatures 
above 10 GK are not reached for the considered trajectory). 

Eor t ^ 1/Atot, Ye reaches the equilibrium value Y^ 


© 2015 RAS, MNRAS 000, [iHT^ 
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given by 

p — 


Ajy + Ag 


A+ + A- 


Aiy + Ag+ + Aj/ + Ag- 


( 8 ) 


The electron (anti)neutrino capture rates can be written in 
terms of the average (anti)neutrino capture cross sections 
(cTj^e/iye) ( Pllunibi et al.|2014 ) as 

Ke> , (9) 


Aj^g — 

Aj?^ — 


47rr2(^j^g) 

Lve 

47rr2(^j^g) 




( 10 ) 


Here, the local (anti)neutrino number densities are ex¬ 
pressed by the ratios of the global luminosities, and 

the mean energies of the radiated neutrinos, mul¬ 

tiplied with the spherical surface 47rr^ that surrounds the 
central neutrino source at a radial distance r (for every 
trajectory we adopt the time-dependent radial distance r 
from our hydrodynamical model). This dilution of the 
neutrino flux is a crude approximation and holds, at best, 
far away from the neutrinosphere, provided the emission is 
isotropic, i.e., if directional variations of the neutrino fluxes 
do not play a role. Close to and below the neutrinosphere, 
however, such a description breaks down but can be justi¬ 
fied by the fact that at these locations electron and positron 
captures dominate and their competition enforces a state of 
weak equilibrium. At large distances the asymptotic electron 
fraction is determined by neutrino and antineutrino absorp¬ 
tions, for which reason direction-dependent differences of the 
neutrino exposure of the ejecta would be important for a de¬ 
tailed discussion of neutrino effects on merger ejecta. Nev¬ 
ertheless, despite these shortcomings, we apply the simple 
ansatz at all radii r where p ^ peq in order to discuss basic 
aspects of the impact of neutrino processes with nucleons in 
the merger ejecta in a parametric way. 

Similar to Eqs. ( IM , the electron and positron capture 
rates are given in terms of the average electron/positron 
capture cross sections (crg-/g+) by 


Ae+ — chg+(crg+), 

Ag- = cng-(crg-). 


( 11 ) 

( 12 ) 


where c is the speed of light and rig- and hg+ the electron 


and positron densities, as detailed in Pllumbi et al. (2014). 

In turn, the average cross sections for electron neutrino 
and antineutrino captures, as well as those for electron and 
positron captures, including the weak magnetism and recoil 
corrections, are taken from Pllumbi et al. (2014) (see also 
[Horowitz Li||1999 ). While the electron and positron cap¬ 
ture rates are temperature- and density-dependent only, the 
(anti)neutrino capture rates depend on the (anti)neutrino 
luminosities and mean energies, hence require a detailed 
knowledge of the neutrino properties at each time step. 

For the present study, we consider representative 
(anti)neutrino luminosities and angle-averaged mean ener¬ 
gies that are assumed to remain constant in time. For a given 
mean energy the electron neutrino temperature 

is deduced from the relation 




Em 

F2(0) 


(13) 


where is the Boltzmann constant and En are the fermi 


integrals (Takahashi et al. 1978) of order n for vanishing 



t [ms] 


Figure 1. (Color online). Time evolution of the electron neutrino 
and antineutrino luminosities from jWanajo et al.j ( |2Q14| ) and the 
corresponding (black solid line) as defined by Eq. 1^. The 

black dotted line gives the without weak magnetism and 

recoil corrections, i.e. = 1 in Eq. The (anti)neutrino 

mean energies are taken consistently from Fig. 1 (lower panel) 
of jWanajo et al.j ( |2014| ). Asymptotic values of Y^^ are calcu¬ 
lated only for non-negligible neutrino luminosities, i.e. for times 
t ^ 4 ms. 



t [ms] 


Figure 2. (Color online). Same as Fig.^ but for the electron neu¬ 
trino and antineutrino luminosities from jRuffert &: JankalpOOTt - 
The (anti)neutrino mean energies are taken consistently from 
Fig. 17 (lower panel; Model Be) of this reference. 


chemical potential, assuming nondegenerate neutrino spec¬ 
tra. A similar expression holds for the antineutrino temper¬ 
ature. For given (anti)neutrino luminosities and mean en¬ 
ergies this allows us to determine all other moments of the 
(anti)neutrino energy spectra, and consequently the corre¬ 
sponding capture cross sections and rates (Eqs. 1^ . 

Assuming that electron (anti)neutrino captures domi¬ 
nate over electron and positron captures and the abundance 
of heavy nuclei is negligible, the asymptotic value of Ye 
(Eq.§ can be approximated by 


yriyoo ^ 


nn 

eJjy, 


T p fmr T fm: 


(14) 
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because the rates can be expressed as oc — 

where = {El)!{Eiy) = F4(0)F2(0)/F3^(0) x {Eiy) 
and corresponds to the weak magnetism and recoil cor¬ 
rections that can be found in [Pllumbi et al. ( |2014| ). The 
corresponding value of is shown in Fig.M for the lumi¬ 
nosities and mean energies of Wanajo et ah |2014 ) and lies 
between 0.25 and 0.50 for times t ^ 5 ms after the merg¬ 
ing of the binary NSs. These values may, however, still be 
modified by effects of heavy nuclei. The weak magnetism 
and recoil corrections on the (anti)neutrino rates are seen 
in Fig. to increase the asymptotic value of by up to 
20% because they reduce the antineutrino capture cross sec¬ 
tion and simultaneously increase the neutrino capture cross 
section ( Horowitz Li||19^ ). 

The neutrino properties, and in particular the antineu¬ 
trino to neutrino luminosity ratio, are found to vary signif¬ 
icantly between different hydrodynamical simulations but 
also depend on the adopted equation of state (Sekiguchi et 


al. 2015). Avoiding the complexity of self-consistent neutrino 
transport in hydrodynamical simulations, we shall restrict 
ourselves in our sensitivity study to constant luminosities 
and mean energies taken at selected times from previous 
simulations ( [Ruffert Janka|2001 Wanajo et al.|2014 ). We 
consider first two representative sets of values for electron 
(anti)neutrino luminosities and mean energies as obtained 
by [Wanajo et al. ( 2014| ) (see their Fig. 1), namely those 

corresponding to the instants of 5 and 6 ms, i.e., 

• Case 1: t ^ 5 ms with = 0.6 x 10^^ erg/s; = 

1.3 X 10^^ erg/s; (E^^) = 12 MeV; {E^,^) = 16 MeV, 

• Case 2: t ~ 6 ms with = 2.6 x 10^^ erg/s; = 

4.0 X 10^^ erg/s; (E^^) = 13 MeV; {E^J = 16 MeV. 

Note that Case 1 leads to an asymptotic value Y^"^ ^ 0.31, 
whereas Case 2 yields Y^"^ ~ 0.42 (Fig. [^. 

In order to test more thoroughly the impact of the neu¬ 
trino processes on the nucleosynthesis, we also consider here 
the electron (anti)neutrino properties calculated in the NS 
merger simulation of jRuffert V Janka] ( 200l| ) (cf Model Be 
in their Fig. 17). The corresponding (anti)neutrino lumi¬ 
nosities are shown in Fig. together with the asymptotic 
values Y^°^ (with and without the weak magnetism and re¬ 
coil corrections). Overall, lower luminosities are found in this 
model in comparison with Wanajo et al. (2014), but also a 


relatively higher emission of electron antineutrinos. Conse¬ 
quently, lower values of Y^^ are predicted for this model. 
In our present sensitivity analysis we also select two cases 
for the electron (anti)neutrino luminosities and mean en¬ 
ergies from the results of Ruffert V Janka] ( |2001| , namely 
those corresponding to times of 6 and 10 ms (and hereafter 
referred to as Cases 3 and 4, respectively), i.e., 

• Case 3 : t ~ 6 ms with = 0.3 x 10 ^^ erg/s; = 

10 ^^ erg/s; {E,^) = 12.5 MeV; {E^,^) = 17.4 MeV , 

• Case 4: t ~ 10 ms with = 0.6 x 10 ^^ erg/s; = 

1.6 X 10 ^^ erg/s; {E^J = 13.5 MeV; = 16.3 MeV . 


Case 3 leads to an asymptotic value Y^"^ ~ 0.21, while Case 
4 yields - 0.29 (Fig. [2]). 

Finally, it should be mentioned that our assumption of 
time-independent (anti)neutrino luminosities can be ques¬ 
tioned, since the neutrino-ejecta interaction is a highly time- 
dependent problem, where the relative time between the 
growth of the neutrino emission and the mass ejection mat¬ 


ters. Considering constant luminosities is a very crude but 
simple approximation, which is sufficiently good to demon¬ 
strate the impact of neutrino processes on the time evolution 
and mass distribution of the electron fraction and the cor¬ 
responding consequences for the r-process. A predictive as¬ 
sessment of neutrino effects on the nucleosynthesis in merger 
ejecta would also have to take account variations of the neu¬ 
trino emission with different directions. Matter expelled to¬ 
wards the polar directions is exposed to different neutrino 
conditions than matter that leaves the system along equato¬ 
rial trajectories. Again, a more detailed description of neu¬ 
trino transport effects is demanded and is beyond the scope 
of our present parametric study. 


3 IMPACT OF /^-INTERACTIONS ON THE 

ELECTRON FRACTION 

3.1 Time evolution of Ye 

To illustrate the impact of electron (anti)neutrino, electron 
and positron captures on the evolution of Ye, we show in 
Figs. and the time evolution for two specific trajec¬ 
tories during the expansion from the initial density peq = 
10^^gcm“^ down to density pnet, where the reaction net¬ 
work calculations are initiated. Both trajectories are studied 
including (anti)neutrino captures with neutrino properties 
for Case 1. 

For trajectory 400720 (Fig. |^, neutrino and antineu¬ 
trino captures on free nucleons with rates up to some 10^ s“^ 
lead to a rapid increase of Ye from 0.17 up to 0.31, corre¬ 
sponding to the asymptotic value of Y^^^ in Case 1. The 
(anti)neutrino rates are found to dominate the electron and 
positron capture rates already at densities slightly below peq- 
This is linked to the slow 1/r^ decrease of the (anti)neutrino 
fluxes, which is more shallow than the steep temperature 
dependence (roughly like T®) of the electron and positron 
capture rates. At t > 4 ms, free neutrons and protons in 
the expanding matter partially recombine into a-particles 
{{Ah) = ^ — 4), thus giving rise to the 

a-effect and therefore a further increase of Ye. The rate A* * 
becomes larger than A+ ~ and approaches {Xh^ + Aj^g)/2 
of material dominated by a-particles. 

For trajectory 486857 (Fig. [^, at the initial density 
Peq, the mass element is characterised by a low temperature 
T < 10^° K and consequently is composed of heavy nuclei 
{Ah) — 120 typical of the low values of Ye — 0.1 in the outer 
NS crust. The fast e~ and (anti)neutrino capture rates lead 
to a rapid increase of Ye , but all these rates are comparable 
and Ye fluctuates wildly. The mass element is then subject to 
a new compression phase that sets in at t ~ 1.5 ms, and the 
corresponding high temperatures (T > 10^^ K) photodisso- 
ciate the matter into free nucleons. It should be noted that 
during this recompression episode, the higher densities cause 
neutrinos to become trapped again, for which reason the 
application of Eqs. ( 9|10 ) remains problematic and highly 
schematic. However, at these high-density, high-temperature 
conditions, electron captures dominate and re-neutronize 
the material until they are counterbalanced by positron cap¬ 
tures, whose rate increases dramatically to achieve weak 
equilibrium with the electron captures. During the subse¬ 
quent expansion phase (at t > 2 ms), Ye rises gradually until 
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Figure 3. (Color online). Upper panel: Time evolution of tem¬ 
perature, density, radius, and mean atomic mass of nuclei heav¬ 
ier than protons, {A^) = ^z>2 ^^1 ^z>2^ ^ between peq and 
Pnet for trajectory 400720. Lower panel: Analogue for the elec¬ 
tron (anti)neutrino, electron and positron capture rates, and for 
A* (Eq. and the electron fraction. The (anti)neutrino proper¬ 
ties correspond to Case 1. The late-time increase of Ye is caused 
by the a-effect, which does not asymptote to a terminal value 
until the network is started. Ye continues to evolve subsequently 
during the nucleosynthesis because of ^-decays. 

(anti)neutrino absorptions take over to push Ye towards its 
asymptotic value of Y^°^ ~ 0.31 (at t ~6ms). With decreas¬ 
ing temperatures, the a-effect finally becomes responsible 
for a late increase of Ye at t ^ 7.5 ms. 


3.2 Ye distributions 

Six different cases are studied to estimate the impact of elec¬ 
tron and positron captures as well as electron (anti)neutrino 
absorption on the ejecta mass distribution as function of 
Ye- In the first case, weak interactions of free nucleons are 
not allowed. In the second case, electron and positron cap¬ 
tures are turned on, but not (anti)neutrino absorptions. In 
the remaining four cases, both electron/positron captures 
as well as (anti)neutrino captures are switched on, with the 
(anti)neutrino properties being defined by Cases 1-4. For 
these six cases, the resulting Ye distributions at density pnet 
are shown in Fig. 

In the first case without weak interactions of free nucle- 


Figure 4. (Color online). Same as Fig. but for trajectory 
486857. 


ons below peq (upper left panel of Fig. [^, the Ye distribu¬ 
tion at pnet is identical to the one given at the initial density 
Peq. This case, however, differs from the standard case we 


considered previously (Goriely et al.||2011 

Bauswein et al. 

2013 Goriely et al. 2013 Just et al.||2015 

) in two aspects: 


First, we start our calculations of weak interactions at den¬ 
sity Peq = 10^^gcm“^ with /3-equilibrium distributions of 
electrons, positrons and neutrinos. This shifts and broadens 
the Ye distribution from previous values of <0.1 to a wider 
range between 0 and ~0.3. Second, no temperature post¬ 
processing is performed here. The initial temperatures of 
the trajectories are significantly higher than those deduced 
from the post-processing applied in our previous studies. 

As visible in Fig. lower left panel, the Ye distribu¬ 
tion at density pnet is significantly affected by -captures 
between peq and pnet, although some low-Ye ( < 0.2) ejecta 
are left. Dominant parts of the ejecta are now found at Y^e 
values between 0.3 and 0.4. When switching on neutrino ab¬ 
sorptions, for any of the Cases 1-4 the asymptotic values 
of Ye^°° are approached, and further enhancement by the 
a-effect produces peaks of the mass distributions in a range 
of Ye values between 0.2 and 0.5. The peak values depend 
sensitively on the adopted neutrino properties. Cases 1 and 
4 lead to rather similar Ye distributions, owing to the fact 
that the neutrino properties are broadly comparable. 

These distributions also depend sensitively on the tem¬ 
perature. When the temperature along the trajectory is ar- 
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Figure 5. (Color online). Histograms of fractional mass distributions of the 4.9 x 10“^Mq of matter ejected in our 1.35-1.35Mq NS-NS 
merger model as functions of Ye at density pnet, assuming that no weak interactions of free nucleons take place below peq (upper left), 
that only electron and positron captures affect the Ye evolution for p < peq (lower left), and including electron and positron captures as 
well as (anti)neutrino captures at p < peq (Case 1: upper centre; Case 2: lower centre; Case 3: upper right; Case 4: lower right). 


tificially increased or decreased by factors of 3, rather dif¬ 
ferent results are obtained, as shown in Fig. using the 
(anti)neutrino properties of Case 1. As before we start the 
expansion evolution at peq with the Ye mass distribution 
shown in the upper left panel of Fig. Reduced temper¬ 
atures diminish the presence of positrons and thus favor 
electron captures compared to positron captures. Moreover, 
lower temperatures also lead to a faster freeze-out of 
captures during the ejection of the mass elements. Without 
neutrino and antineutrino absorptions, lower temperatures 
therefore tend to neutronize the ejecta for p < peq and the 
mass distribution becomes more narrow and is shifted to 
lower values in the range of 0 < Fe ^0.1 (Fig.j^ upper left 
panel). Including neutrino and antineutrino absorptions, re¬ 
duced temperatures have the opposite effect in pushing the 
mass distribution to higher values of Ye (with a peak above 
0.4) compared to the standard-temperature result for Case 1 
in Fig. (upper right panel, with a peak of the distribution 
between Ye — 0.2 and 0.3). This behavior can be understood 
by the efficient recombination of free nucleons to a parti¬ 
cles and heavy nuclei, which strengthens the heavy-nuclei 
(a) effect so that the peak of the distribution wanders to 
Ye — 0.45. On the other side, increased temperatures re¬ 
duce the electron degeneracy and thus allow for the pres¬ 
ence of higher positron densities, thus enhancing positron 
captures on neutrons. In addition, captures continue for 


a longer period of time along the ejecta trajectories. With¬ 
out (anti)neutrino absorption, these effects shift the Ye mass 
distribution from the initial one at peq (upper left panel of 
Fig.[5f towards higher values of Ye. This shift is stronger for 
more slowly expanding mass elements and weaker when the 
expansion is very fast. Correspondingly, the mass distribu¬ 
tion versus Ye at pnet is very broad and stretches from ~0.03 
up to a very pronounced maximum close to 0.5, because 
the capture equilibrium at high-entropy conditions favors 
symmetric conditions with respect to neutrons and protons. 
Taking into account (anti)neutrino absorption prevents this 
dramatic shift towards Ye ~ 0.5, because at large distances 
neutrino captures dominate absorptions and therefore Ye 
asymptotes to values around Y^"^ (~ 0.31 for Case 1, lower 
right panel of Fig.[^. Since the high temperatures favor nu¬ 
cleons and suppress the early formation of a particles and 
heavier nuclei, the influence of the a effect is clearly weaker 
than in the case of reduced temperatures (compare lower and 
upper right panels of Fig.[^. As mentioned above, our calcu¬ 
lations with varied temperatures are not consistent with the 
initial Ye distributions used at a density of peq, because these 
distributions are calculated for the original ejecta tempera¬ 
tures provided by the hydrodynamic NS-NS merger model. 
However, Fig. [^demonstrates that asymptotic values deter¬ 
mined by neutrino capture equilibrium and influenced by 
the a (heavy-nuclei) effect are reached for most of the tra- 
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Figure 6. (Color online). Histograms of fractional mass distribu¬ 
tions as functions of Ye at density pnet for captures exclusively 
at p < peq (left) and including electron, positron captures as well 
as (anti)neutrino captures during the evolution from peq to pnet 
(Case 1; right), when the temperatures are decreased (upper pan¬ 
els) or increased (lower panels) artificially by factors of 3. 


jectories. The final Ye values can therefore be expected to 
mostly have lost the memory of the initial conditions at peq. 
Correspondingly, when neutrino absorptions are included, 
the final mass distributions of Ye (at pnet) are considerably 
more narrow than the relatively broad distribution of initial 
Ye values before the expansion from peq to pnet (compare the 
right panels of Fig. with the upper left panel of Fig. |^. 


4 R-PROCESS NUCLEOSYNTHESIS 


At densities p < pnet, the abundance evolution is followed 
by a full reaction network (for more details, see Goriely et 


al.|2011||Bauswein et al.| 2013)[Goriely et al.|2013||Just et al. 

2015| , which in contrast to earlier studies includes now also 
the weak interactions of free nucleons as detailed in Sect.|2] 
The final abundance distributions are shown in Fig.j^for the 
same six different cases discussed in Sect. |3.2] and displayed 
in Fig. namely for a case without weak interactions of 
free nucleons below peq, a case including -captures but 
without (anti)neutrino absorption reactions, and four cases 
where both e^- and (anti)neutrino captures are taken into 
account at p < peq (Gases 1-4 for the v properties). 

Because of its low-Ye distribution, the case neglecting 
/3-interactions gives rise to an elemental abundance distri¬ 
bution similar to the one obtained in our previous studies 


(Goriely et al.|20lT Bauswein et al. 2013 Goriely et al. 2013 


Just et al.|2015 ). It is characterized by the production of es¬ 
sentially only A > 140 nuclei through several loops of fission 
recycling. When -captures on free nucleons are switched 
on, some low-Ye ( ^ 0.2) material can still lead to nucleosyn¬ 
thesis with fission recycling and a significant production of 
the third r-process peak, but the higher-Ye (0.3-0.4) matter 
can now also contribute to the production of 90 ^ A ^ 140 
nuclei with a strong second iV = 82 peak. 

When (anti)neutrino captures are included, too, nucle¬ 
osynthesis with fission recycling does not take place any 
longer, but the final abundance distribution still resembles 


the one in the solar system fairly well. A significant amount 
of 50 ^ A ^ 90 nuclei can now also be produced in addition 
to the A > 90 r-nuclei, especially for (anti)neutrino prop¬ 
erties corresponding to Gase 2, as shown in Fig. In this 
case, important element production around ®°Ni is obtained, 
originating from ejected mass elements with Ye ^ 0.45. 

Without yd-interactions below peq, about 90% of the 
ejected material is found to be r-process rich. If -captures 
on nucleons are switched on, the material is made of 76% 
r-process nuclei, and when (anti)neutrino captures are effec¬ 
tive, too, we find that 45% of the ejected matter is made of 
r-nuclei in Gase 1, only 1.9% in Gase 2, 67% in Gase 3 and 
40% in Gase 4. In the last four cases, a significant part of the 
material is made of a-particles (~ 22% in Gases 1 and 4, 50% 
in Gase 2, and 17% in Gase 3), and the remaining part con¬ 
sists of 50 ^ A ^ 70 nuclei (see Fig. [^. Despite the rather 
robust production of a solar-like distribution of r-nuclei, the 
absolute and relative amounts of r-material vary strongly 
from case to case, sensitively depending on the neutrino ex¬ 
posure of the ejecta as determined by the assumed properties 
of the neutrino emission. Neutrino properties corresponding 
to Gase 2 yield a total amount of ejected r-material that is 
significantly smaller than for the other three sets of neutrino 
properties. 

The total, mass-averaged nuclear energy-release rate 
that is available for heating the ejecta per unit of mass, the 
average temperature of the ejecta, and the average atomic 
mass number of the ejected abundance yields are plotted as 
functions of time for our six studied cases in Fig. After 
some 10 s and up to nearly one day, the energy release rates 
of all cases are fairly similar except for the treatment of neu¬ 
trinos and /3-interactions according to Gase 2. In this case, 
the large electron fractions of the ejecta favor the produc¬ 
tion of light elements and only a small amount (~ 2%) of 
the ejecta material consists of r-process nuclei (Fig. and 
is therefore subject to fast ^d-decays. At late times, typically 
after one day, an additional source of energy is found from 
the a-decay of long-lived heavy nuclei that can only be sig¬ 
nificantly produced when no yd-interaction of nucleons are 
included below peq. Despite the considerable differences of 
the nuclear energy-release rates at early times, the cooling 
evolution as measured by the mass-averaged temperatures 
is rather similar in all cases. 

The time evolution of the average nuclear mass num¬ 
ber of the ejected material clearly shows that only without 
/3-interactions of free nucleons a significant amount of fissile 
nuclei can be produced with (A) reaching values up to 170, 
which is slightly below the value of about 200 obtained when 
the Ye distribution of the inner crust of a cold NS is consid¬ 
ered as initial state (see, in particular, [Goriely et al.||2011 ). 
The sequence of decreasing values of (A) follows basically 
the hierarchy of increasing values of the mean electron frac¬ 
tion, namely (Ye) =0.11 when /3-interactions are ignored, 
(Ye) =0.27 when only -captures are taken into account, 
and (Ye) = 0.25, 0.34, 0. 35, 0.45 in the Gases 3, 1, 4 and 2, 
respectively. The only (slight) inversion is obtained for the 
calculation with only e^-captures, because in this case the 
Ye distribution is very wide (Fig. lower left panel) and 
the significant amounts of low-Ye (~ 0.05 — 0.15) material 
contribute to the production of heavy nuclei which increase 
the average nuclear mass number (A). 
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Figure 7. (Color online). Abundance distributions of the 4.9 X 10“^Mq of material ejected in our 1.35-1.35Mq NS-NS merger model, 
assuming that no weak interactions of free nucleons take place at p < peq (upper left), only electron and positron captures happen below 
Peq (upper right), and electron and positron captures as well as (anti)neutrino captures affect the Ye evolution for p < peq (Case 1: 
upper centre; Case 2: lower centre; Case 3: upper right; Case 4: lower right). The abundance distributions are normalized to the solar 
r-abundance distribution (open circles) in the rare-earth peak {A = 165) region. 



Figure 8. (Color online). Mass fraction as a function of 
the atomic mass for the matter ejected by our 1.35-1.35Mq 
NS merger model, including electron, positron as well as 
(anti)neutrino captures for Cases 1 and 2 of the neutrino proper¬ 
ties. 

5 SUMMARY AND CONCLUSIONS 

In this paper we reported the results of a parametric study 
to investigate how /3-interactions of free nucleons can affect 
the Ye evolution and mass distribution in NS merger ejecta 



Figure 9. (Color online). Time evolution of the total radioactive 
heating rate per unit mass, {Q) (red), mass number {A) (blue), 
and temperature (T) (black), all mass-averaged over the ejecta, 
for the 1.35-1.35Mq NS merger. The solid lines correspond to the 
case without /3-interactions of free nucleons at p < peq, the dotted 
lines to the case where only -captures are taken into account, 
the dashed lines to the case where and (anti)neutrino captures 
for neutrino properties according to Case 1 are included, the dash- 
dot lines for neutrinos according to Case 2, the long-dash-dot lines 
correspond to Case 3, and the double-dash-dot lines to Case 4. 
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and the corresponding nucleosynthesis. To this end we used 
the temperature-density trajectories of a large set of mass 
elements representing the ~ 5 x 10“^ M© of matter ejected 
in a relativistic merger simulation of a symmetric l.SSM©- 
l.SSM© NS binary with the non-zero temperature DD2 nu¬ 
clear equation of state. Using the total lepton number pro¬ 
vided by our NS-NS merger simulation, we assume matter to 
be in /3 equilibrium at the temperature of the hydrodynami- 
cal model and a fiducial density of peq = 10^^ gcm“^. These 
conditions define the starting points of our post-processing 
of the composition histories of the considered ejecta ele¬ 
ments, for which we take into account electron and positron 
captures as well as electron neutrino and antineutrino ab¬ 
sorption on free neutrons and protons. Below a density pnet, 
where the temperature has decreased to 10^° K (or the neu¬ 
tron drip density, if matter below peq remains cooler than 
10^° K) the full network calculation is applied instead of nu¬ 
clear statistical equilibrium. Our description of weak inter¬ 
actions of free nucleons includes weak magnetism and recoil 


corrections according to Horowitz & Li (1999) and Pllumbi 


et al. (2014). Avoiding the complications of treating neu¬ 


trino transport, we simply use exemplary data from publi¬ 
cations for prescribing the neutrino luminosities and mean 
energies needed to compute the neutrino absorption rates 
in our parametric approach (“Cases 1-4”). This elementary 
prescription also accounts for the still large uncertainties of 
the model predictions for the neutrino emission and its di¬ 
rectional asymmetries in the generically three-dimensional 
merger scenario. 

Our modeling strategy follows the spirit of previous, nu¬ 
merous parametric investigations of nucleosynthesis in the 
neutrino-driven wind of newly formed neutron stars in su- 


lArnould et al.| 
therein) and of; 
of compact ob, 

|20071 IPllumbi et al.||2014[ and references 
accretion tori around black holes as remnants 

ject mergers (e.g., Surman h McLaughlin 

2005 

ISurman et al.||2008| | Wanajo & Janka|2012| |Caballero 


et al.||2012). In particular, we track in detail the charged- 


current /3-interactions of neutrinos with free nucleons, which 
determine the evolution of the electron fraction outside of 
the neutrino-trapping regime, and employ a full set of trajec¬ 
tories that characterizes the conditions in dynamical ejecta 
from the merging phase of a representative binary neutron 
star. These conditions differ from proto-neutron star and 
accretion-torus winds not only concerning the range of en¬ 
tropies. The ejecta, especially, possess much faster expansion 
velocities, which for the bulk of the matter can be 25-50% of 
the speed of light, for some fraction of the ejecta even faster, 
whereas neutrino-driven proto-neutron star winds have typi¬ 


cal velocities of 3-7% of the speed of light (e.g., Arcones et al. 


2007), neutrino-driven winds from massive neutron stars as 


relics of NS mergers may achieve expansion velocities up to 


8-10% of the speed of light (Perego et al. 2014), and the main 
mass of neutrino-driven outflows from BH-accretion tori can 
reach 10-20% of the speed of light ( Just et al.||2015 ). The 
correspondingly shorter expansion time scales of dynamical 
merger ejecta can enable a strong r-process even for moder¬ 
ately low Ye- 

Our calculations confirm recent results of [Wanajo 
al. (2014) that solar-like r-process abundances are produced 
in dynamical NS merger ejecta even when ^d-reactions of 
free nucleons are taken into account and lead to a signif¬ 


icant increase of the average electron fraction. In contrast 
to previous works, where charged-current neutrino-nucleon 
interactions were ignored, however, also nuclei with mass 
numbers A < 140 are ejected in larger amounts. 

In detail, our results can be summarized by the follow¬ 
ing points: 


• Ignoring /3-interactions at densities p < peq, we confirm 


our previous results of 

Goriely et al. (2011); Bauswein et al. 

o 

CO 

Goriely et al. ( 

2013); Just et al. (2015) that almost 


exclusively r-nuclei in the regime A > 140 are produced, in 
spite of a moderate increase of the average ejecta Ye asso¬ 
ciated with the assumption of /3-equilibrium at density peq 
instead of our previous use of electron fractions of cold neu¬ 
tron star crust matter. 

• Positron and electron-neutrino captures, enhanced by 
weak magnetism corrections (which increase the absorption 
cross section of i^e and reduce that of Pe) and supported 
by the a effect, lead to a shift of the average ejecta Ye to¬ 
wards higher values when matter expands downwards from 
an initial density of peq = 10^^ gcm“^. 

• Captures of cause a wide spread of the Ye mass dis¬ 
tribution, reaching from values of Ye 0.1 up to 0.4-0.5. 
This reflects the wide range of thermodynamic conditions of 
the ejecta at peq with ejecta trajectories that describe cool 
as well as hot conditions. Absorption processes of i^e and 

take over when the temperature of the expanding mat¬ 
ter has fallen to low values where absorptions become 
ineffective. These (anti)neutrino captures tend to push Ye 
towards the asymptotic value YJ^°° for i^e-i^e capture equi¬ 
librium. The corresponding Ye mass distributions (at pnet, 
where the full network calculation was started) are rather 
narrow in all cases, with a spread of Ye values of only 0.1- 
0.15. The mean values of the distributions, however, vary 
between ~0.25 and ~0.45, depending on the relative size of 
the i^e and z^e luminosities. The latter are highly uncertain 
and variable and are sensitive to the binary-parameter de¬ 
pendent merger dynamics, the nuclear equation of state, and 
the still not well determined directional asymmetries of the 
neutrino emission. 

• In the presence of yd-interactions of free nucleons the 
temperature also plays an important role for the Ye evolu¬ 
tion of the ejecta. It is another aspect of the hydrodynam- 
ical merger models that depends on the system properties 
and the detailed ejection dynamics (which differ between 
different ejecta components) and can be numerically prob¬ 
lematic, because resolution and numerical/artificial viscosity 
can have an influence on the accuracy of the determination 
of the thermal conditions. Lower temperatures reduce the 
effects of captures but increase the importance of the ot 
effect (thus pushing the average Ye closer to 0.5), whereas 
higher temperatures enhance captures but nevertheless 
have little influence on the final Ye distribution at pnet, for 
which the asymptotic value of Ve-Ve capture equilibrium is 
more relevant than the initially fast captures. 

• In all investigated model cases, the production of heavy 
r-process matter with strong second and third abundance 
peaks and a near-solar distribution in the rare-earth region 
is a robust outcome. However, in contrast to previous results 
where weak interactions of free nucleons were ignored, also 
considerable amounts of matter are synthesized to A < 140 
nuclei. The strength of the production of A < 90-100 ma- 
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terial is sensitive to the neutrino-emission properties that 
determine the (anti)neutrino absorption. In extreme cases 
where Ye gets close to 0.5, significant amounts of iron-group 
nuclei (A ~ 50-60) can be ejected. The relative fraction of 
heavy r-process matter (from the second peak upward) in 
the ejecta therefore varies dramatically between the differ¬ 
ent investigated cases of neutrino-emission conditions and 
spans a range from more than 75% down to just 2%. 


We emphasize that our parametric approach is highly 
simplified and ignores important neutrino-transport effects 
like the exact spectral distribution of the neutrino fluxes, 
the direction dependence of the neutrino emission and corre¬ 
sponding precise radial dilution function, and the time evo¬ 
lution of the neutrino emission relative to the ejection time 
of the matter. Nevertheless, our results, in support of those 


of Wanajo et al. (2014) and Sekiguchi et al. (2015), have 
important consequences for the further exploration of the 
nucleosynthesis connected to compact binary mergers and 
the discussion of astrophysical implications. In fact, they 
demand a major revision of the current picture of r-process 
production in such events. 

In view of our results it is obvious that a proper treat¬ 
ment of the neutrino physics, in particular of the neutrino 
irradiation of the ejected material, is essential for making 
quantitative predictions of the elemental yields and espe¬ 
cially of the total mass of r-process material that is thrown 
out by the dynamical merger ejecta. Approximations like the 
ones used in our study can be satisfactorily removed only 
when ultimately detailed, three-dimensional neutrino trans¬ 
port is consistently included in the hydrodynamical simula¬ 
tions. 

Our study suggests that the relative contributions of 
matter with A < 90, 90 < A < 140 and A > 140 are likely 
to depend strongly on the binary properties and even the di¬ 
rection of mass ejection, in addition to the equation of state 
dependence that was found for the electron-fraction distri¬ 


bution in the recent work of Sekiguchi et al. (2015). Differ¬ 


ent from expectations so far, this means that symmetric or 
nearly symmetric binary NS mergers could exhibit a signif¬ 
icantly different ejecta composition than highly asymmetric 
NS-NS mergers and NS-BH mergers, in which the NS is dis¬ 
rupted before it can be swallowed by the BH. In the last two 
cases the lower-mass component develops an extended tidal 
tail, from which considerable amounts of cold, unshocked 
matter can be centrifugally ejected before the neutrino lu¬ 
minosities rise high and thus before neutrino exposure of 
these ejecta plays an important role. In such a situation 
the ejecta will not only be expelled highly anisotropically 
but will also carry a far dominant fraction of the mass in 
the form of A ^ 140 nuclei as predicted in previous studies 


(e.g., Goriely et al. poTT | [Bauswein et al.|2013||Goriely et ah 


[2013 Just et al.||2Q15 ). In contrast, in symmetric or nearly 

symmetric NS mergers the contribution of A ^ 140 material 
will be higher. If the merger remnant collapses to a BH on a 
millisecond time scale, neutrino exposure of the ejecta may 
be avoided, but the broad Ye distribution caused by cap¬ 
tures will allow for the production of A > 90 nuclei with a 
strong second A" = 82 peak. If the merger remnant remains 
transiently or permanently stable, neutrino exposure of the 
dynamically expelled matter becomes important, enabling a 
higher production of A < 90 species. For extreme cases of a 


luminous Ue flux, the ejecta might then even be dominated 
by iron-group nuclei including radioactive nickel isotopes. In 
particular, however, the exact composition and the relative 
fraction of high-mass and low-mass species could depend on 
the direction of the mass ejection. If most of the neutrino 
flux is emitted to the polar directions due to the rotational 
deformation of the merger remnant, matter expelled near the 
equatorial plane will receive less neutrino exposure in addi¬ 
tion to its potentially faster escape. This will allow for more 
neutron-rich conditions close to the equator whereas polar 
ejecta may contain more proton-rich contributions. Future, 
more complete merger models with neutrino transport will 
have to clarify these possibilities. 

Since the photon opacity, k, of the expanding gaseous 
ejecta is strongly dependent on the presence of high-opacity. 


et al.|2013 Tanaka & Hotokezaka|2013 Tanaka et al. 

2014), 

the relative contribution of trans-iron eleme 
will have a severe impact on the peak lum 
peak time, tpeak oc and the effe 

perature, Tpeak oc of the electroma 

that is expected from the radioactively he 
ejecta cloud (“macronova” or “kilonova”. 

nts to the ejecta 
inosity, Tpeak oc 
active peak tem- 
.gnetic transient 
ated, dynamical 
Li & Paczyhski 

1998 


|Kulkarni||2005| |Metzger et al.||2010 

IRoberts et al. 


2011 


Goriely et al. 112011). The current picture of merger- 


type and remnant-type dependent redder (near-infrared) or 
bluer emission ( Metzger &: Fernandez|[^14[ [Perego et al. 
2014) and in particular of the envisioned late-time infra-red 


radiation component from the dynamical ejecta (Kasen et 
aL]|2014 ) might require revision or extension in view of our 
results. 

Finally, it is evident that the strong impact of neutrinos 
on the neutron-to-proton ratio and the nuclear composition 
of the dynamic merger ejecta must move neutrino oscilla¬ 
tions into the focus of interest. It will be necessary to study 
the effects of collective neutrino oscillations (see, e.g., Duan 
et al.||2010 for a review) with similiar intensity as this sub¬ 


ject currently receives in the context of the neutrino emission 
from newly born neutron stars in supernovae. 
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